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ABSTRACT 

Feasibility studies for the conversion of the Nigerian Research Reactor- 1 (NIRR-1) have been 
performed using WIMS and CITATION codes at the Center for Energy Research and Training (CERT), 
Ahmadu Bello University, Zaria Kaduna State. In this work, the core neutronics calculation of NIRR-1 
concerning mass loading of U-235 in the core, shut down margin (SDM), safety reactivity factor (SRF), 
control rod worth, and control rod critical depth of insertion were investigated at low enrichment. Two 
fuel types (UA1 4 and UO2) were considered and the uranium densities required for the conversion of 
NIRR-1 core to low enrichment were computed to be 1201g/cc with 20% enrichment, 1144g/cc with 
19.75% enrichment, 1274g/cc with 15% enrichment, 1448g/cc with 10% enrichment for UA1 4 fuel type 
and 1141g/cc with 20% enrichment, 1144g/cc with 19.75% enrichment, 1216g/cc with 15% 
enrichment, and 1389g/cc with 10% enrichment for U0 2 fuel type. Significantly, higher uranium 
densities are required to convert NIRR-1 from HEU to LEU - indicating a drastic review of the NIRR-1 
core. 
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INTRODUCTION 

As part of the on-going global effort to convert HEU reactor cores to Low Enriched Uranium (LEU) cores under 
the Reduced Enrichment for Research and Test Reactors (RERTR) program, there is need to study the 
possibilities of converting NIRR-1 core to less than 20% enrichment. Works have been done on Research and 
Test Reactor core conversion around the globe. In some proposed models, the total number of fuel pins (Khamis 
and Khattab, 1999) and the core radius/height ratio (Matos and Lell, 2005) has been drastically changed. This 
brings about noticeable changes in the relative flux values for both inner and outer irradiation sites. In this work, 
HEU and LEU cores are analyzed using the present UA1 4 fuel and a potential LEU fuel (U0 2 clad in zircalloy) 
are considered. The existing HEU core was also analyzed to validate the reactor model used. A significant 
feature of this work is the preservation of the technical and the geometric specification of the reactor so as to 
maintain the original designed thermal hydraulic of the reactor. A detailed description of the HEU core of 
NIRR-1 can be found in the Final Safety Analysis Report (FSAR); (Azande, et al, 2010; SAR, 2005). 

THEORITICAL CONSIDERATION 

The resonance treatment in WIMS is based on equivalence theorems, which relate a library of resonance 
integrals for each resonance absorber in each group to a particular heterogeneous problem. The equivalence 
theorems employed take account of non-homogeneous moderators which might be mixed with the fuel such as 
oxygen or carbon; (Askew et al, 1966). A first order correction for the interaction associated with the presence 
of several resonant nuclides was also included. 

The library of homogeneous resonance integrals has been formed through the use of the Spin Dependence 
Recombination (SDR) code; (Brissenden and Durston, 1965). The SDR code solves by numerical methods the 
slowing down problem associated with moderators homogeneously mixed with resonant nuclides. The code 
solves the Chernick and Vernon integral equation (Rothenstein, 1960) in the form 
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P 0Q {E) is the self collision probability in fuel region, F„{E") is the collision density, P^E) is collision 
density for a general moderator region /, and P ic . is the probability of a neutron scattered in region i suffering its 
next collision in region o. 

In order to derive the equivalence principle, it is usual to assume that the resonances are narrow and well 
separated and the fluxes fp. (fj ]) for energies outside the resonances may be represented by 

*,(«} = p (2) 

In which d^ is a disadvantage factor, which is assumed invariant with energy such that 

If it is further assumed that the probabilities P it> {Ej can be represented by the type of approximation attributed 
to (Bonalumi, 1961) in the form 

zj =-£—£. a = f-i _ p 1« _£__£ (•#* 

to jfT 'ci v -1- 'fp/^cii^v V ,J 

C p j is the probability of those neutrons escaping collision in region o suffer their next collision in region i. 
Applying equations (2), (3), and (4) in equation (1) with the narrow resonance approximation gives 



$M = 



f.^ + Ci-oJ^G* 






I w 



In the case where d^is unity and 2i; G p j — I, gives 

AE 



J = J[^^+(l-0 



*«*-£ (6) 



Z, op is potential scattering cross-section of the fuel with admixed moderator and er^ is the resonant 

macroscopic absorption cross-section. 

For isolated bodies, the Bell modification (Bell, 1959) of the Wigner relational approximation gives 

f Z + oS dE 

a is the bell factor. The methods adopted or extending the equivalence theorem for isolated rods to infinite 
arrays were suggested by Leslie and Jonsson, 1964. 

METHODOLOGY 

A simplified scheme of WIMS and CITATION codes (Azande et al, 2010; Balogun, 2003) were employed to 
generate group constants for the fuel region and non fuel regions such as the control rod, control rod follower 
region as well as reflectors. Unit cell calculations based on Wigner-Seitz cell modeling was used. The core 
lattice was represented in the form of a super cell; (Iqbal et al, (2002). A two-dimensional CITATION input 
(that is r-z cylindrical geometry) (Balogun, 2003) was prepared with group constants associated with 
homogenized regions, including the control rod. Although the group constants associated with the control rod 
was incorporated, the base input data was prepared such that water completely occupies the channel in which 
the rod travels. Data were introduced on four levels including 
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1. Reactor level: data used to describe general calculation options such as 

• the number of components forming the reactor, 

• the number of energy groups used in the calculations, 

• the total number of materials of which the reactor components are formed, 

• some factors indicating to the system how to model the reactor angularly. 

2. Component level: data used to describe the component itself such as 

• the component number; 

• the number of pieces forming the component; 

• for each piece forming the component data include the following parameters: 

- the internal radius of the piece (mm), 

- the external radius of the piece (mm), 

- the height of the piece (mm). 

- some other parameters that can represent holes and other forms of the piece other than the 
cylindrical one. 

3. Material level: data used to describe the elements forming the material. This level is mainly directed to 
WIMSD4 to enable the generation of the group constants for the different materials. 

4. CITATION level: the data here are free format numbers representing the user options to be applied in 
CITATION code during core calculations. This level includes also control cards options of the code 
CITATION. CITATION input file would automatically be compiled based on the indications of the data 
presented here. 

Once the reactor pieces, materials, and CITATION options are specified, an automatic processing of these data 
is performed, and a modeled reactor is generated. When the modeled reactor is generated, material data and 
CITATION options are used to generate cross sections and group constants by the automatic running of 
WIMSD4. Once WIMS is run, an automatic compilation of the CITATION input file is performed; hence 
CITATION runs automatically to generate an output file with the options specified in the CITATION file. 

The initial excess reactivity p g:e is computed automatically as 

and the result is displayed on the screen, ft , f is the effective multiplication factor of the reactor. Flux and 

power distributions are also automatically generated. The generation and transfer of the group constants form 
WIMSD4 input file in WIMS format to the input file of CITATION and in CITATION format is achieved 
automatically. The model used in WIMSD4 to generate the Group Constant is automatically selected too. 

RESULTS AND DISCUSSION 

The results from Tables 1 and 2 show SDM values within the range of 2.0 and 3.5 mk which is in good 
agreement with design specification for shutdown margin of MNSR. The excess reactivity values obtained are 
also in good agreement with the design specifications of (3.5 - 4.0) mk; (Chengzhen, 1993). 
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Table 1: Results for UAL fuel at low enrichments 



Enrichment 


235 T t 

U core 


S.R.F 


Rod worth 


S.D.M 


Excess 


Critical Depth 


(%) 


loading (g) 




(mk) 


(mk) 


Reactivity (mk) 


(cm) 


05.00 


2505 


1.24 


4.36 


0.84 


3.52 


18.69 


10.00 


1448 


1.59 


5.63 


2.08 


3.56 


13.66 


15.00 


1274 


1.70 


5.99 


2.46 


3.53 


12.94 


20.00 


1201 


1.60 


6.16 


2.18 


3.86 


13.66 



Table 2: Results for UOi fuel at low enrichments 



Enrichment 


235 U core 


S.R.F 


Rod Worth 


S.D.M 


Excess 


Critical Depth 


(%) 


loading (g) 




(mk) 


(mk) 


Reactivity 
(mk) 


(mm) 


05.00 


2505 


1.24 


4.36 


0.84 


3.52 


18.69 


10.00 


1389 


1.67 


5.87 


2.35 


3.52 


12.94 


15.00 


1216 


1.77 


6.27 


2.73 


3.55 


12.22 


20.00 


1141 


1.83 


6.48 


2.93 


3.55 


12.22 



When the value of multiplication factor k B ff is greater than unity, the reactor assembly is said to be 

supercritical and when less than unity, it is a subcritical assembly. For fission chain reaction to be just 
continuing at a steady state, k g ff has to be unity. Expectations are that when the control rod is fully withdrawn, 

the super-criticality should be attained while sub-criticality should be attained when the control rod is fully 
inserted. It is therefore, desirable that criticality should be attained about midway between full insertion and full 
withdrawal of the control rod so as to allow for easy control and regulation of reactivity. The results obtained 
show critical depths which approximately satisfied this expectation. The safety reactivity factor (SRF) of greater 
than 1.5 was obtained which is in agreement with Balogun, (2003). 

CONCLUSION 

In this core neutronics calculation of NIRR-1 at LEU core using WIMS and CITATION codes, results are in 
good agreement with HEU results and design specifications. 

ACKNOWLEDGEMENT 

The author acknowledges the management and staff of Center for Energy Research and Training, the Head and 
staff of Physics Department, Ahmadu Bello University Zaria for their various contributions towards the success 
of this work. 

REFERENCES 

Askew J.R, Fayer F.J. and Kemshell P.B. (1966). A General Description of Lattice Code WIMSD Journal of the 

British Nuclear Energy Society. 5: 564 

Azande, T. S. Balogun, G.I. Ajuji, A.S. Jonah, S.A., Ahmed, Y.A. (2010). The use of WIMS and citation codes 
in fuel loading required for the conversion of HEU MNSR core to LEU. Annals of Nuclear Energy , 1223-1228. 

Balogun G.I. (2003). Automating some analysis and design calculations of miniature neutron source reactor at 
CERT (I). Annals of Nuclear Energy 30: 81-92. 

Bell, G.I. (1959). Theory of effective cross-sections. LA-2322, 

Bonalumi, R. (1961). Neutron first collision probabilities in reactor physics, Energia Nuci, Milano, 8, 326. 

Brissenden, R.J. and Durston, C. (1965). The calculation of neutron spectra in Doppler region. Conference on 
the Application of Computing Methods to Reactor Problems, ANL, 7050:51 

Chengzhen, (1993). Experiment adding top beryllium shims for MNSR. CIAE Technical Report. 

Halsall M.J. (1980). A summary ofWIMSD4 input options. AEEW-M137, UKAEA. 



45 



Azande T. S: Continental J. Applied Sciences 6 (1): 42 - 46, 201 1 



Iqbal, M., Abdullah, M., Pervez, S. (2002). Parameteric tests and measurements after shimming of a beryllium 
reflector in a miniature neutron source reactor (MNSR) Ann. Nucl. Energy 29, 1609-1624. 

Khamis, I. and Khattab, K. (1999). Lowaring the enrichment of the Syrian miniature neutron 
source reactor, Annals of Nuclear Energy, 26: 1031-1036. 

Leslie, D.C. and Jonsson, A. (1964). Improvements to the theory of resonance escape in heterogeneous fuel.l. 
Regular Arrays of fuel Rods. AEEW-R 353. 

Matos, J.E., Lell, R.M., (2005). Feasibility study of potential LEU fuels for a generic MNSR reactor. In: 
Proceedings of the 27th International Meeting on Reduced Enrichment for Research and Test Reactors 
(RERTR), Boston, USA, November 6-10. 

Rothenstein, W. (1960) Collision probabilities and resonant integrals for lattices, Nucl. Sci. Engng, 7:162 

SAR, (2005). Final safety analysis report of Nigeria Research Reactor- 1, CERT Technical Report-CERT/NIRR- 
l/FSAR-01. 

Received for Publication: 04/02/10 
Accepted for Publication: 28/04/1 1 



46 



